Critical Current Peaks at 3£>$ in Superconductors with Columnar Defects: 

Recrystalizing the Interstitial Glass 
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The role of commensurability and the interplay of correlated disorder and interactions on vortex 
dynamics in the presence of columnar pins is studied via molecular dynamics simulations. Simula- 
tions of dynamics reveal substantial caging effects and a non-monotonic dependence of the critical 
current with enhancements near integer values of the matching field B$ and SB^ in agreement with 
experiments on the cuprates. We find qualitative differences in the phase diagram for small and 
large values of the matching field. 

PACS numbers: 74.25.Qt, 75.40.Mg, 67.40.Yv, 67.40.-w 
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The rich phases of vortex matter in high temperature 
superconductors result from a complex interplay between 
disorder, interactions and fluctuations^]. Understanding 
the physics of these competing effects has been a chal- 
lenge which underlies the nature of various static and 
dynamic phases and transitions between them. Of par- 
ticular interest is the role of columnar defects created by 
heavy ion radiation in forming more effectively pinned 
phases of interacting vortices. Extended or correlated 
disorder is much more effective at pinning vortices than 
uncorrelated (point-like) disorder and produces upward 
shifts of several tesla in the irreversibility lines and con- 
comitantly larger critical currents. 

The physics of 3D vortex pinning via correlated disor- 
der can be approximately mapped onto the problem of 
interacting quantum bosons in the presence of uncorre- 
lated disorder in 2D 2j. This has formed the basis for 
many computational studies via quantum Monte Carlo 
simulations of local repulsive bosons Q and exact diag- 
onalization studies^]. Several pinned phases of vortex 
matter emerge as a function of B/B&, where J3$ is the 
equivalent matching field where the number of vortices 
equals the number of defects. For the case where the pins 
outnumber the vortices (B < -B$), a Bose-glass phase is 
formed where vortices are localized onto columnar de- 
fects and possess an infinite tilt modulus for tilts away 
from columnar alignment. For equal numbers of vor- 
tices and defects (B = £>$), an analog to a Mott insu- 
lating phase exists which possesses an infinite compres- 
sion modulus 0- The existence of the Bose-glass phase is 
well documented 0, and evidence for the Mott phase 
has been found in recent magnetization relaxation mea- 
surements at low temperatures |5|, |6(. For the case when 
the vortices outnumber the pins, (B > £?$), the sit- 
uation is much less clear. Recent work has suggested 
that vortices not accommodated to the columnar defects 
are caged by pinned vortices to form a weakly pinned 
"interstitial Bose" glass. It has been conjectured that 
the melting temperature T m decreases rapidly as the 
density of interstitial vortices increases for large fields, 
and either extends smoothly into the Bose glass phase 



or shows a change in slope (kink) at B = 0|. 
Some experiments 0, 0, 0, do not show substan- 
tial changes of the irreversibility line at B ^ , w hile others 
show only a mild kink at d H [H Hi The 
change is even less significant for larger defect concentra- 
tions (larger B$). 

Moreover, interesting effects are seen as B further in- 
creases past B$. In analogy with the presence of Mott 
lobes in the dual superconducting-insulator transition in 
3D, one would expect different phases at commensurate 
values of vortex filling. Nowak et al found that the crit- 
ical current in Tl2Ba2Cu06+<5 decreases for B > B^ but 
then increases and goes through a maximum near 3-B^ ja| , 
which is currently unaddressed by theory. 

One of the weakness of previous numerical studies |j 
is the use of a short-range screened interaction which 
misses the softening of the shear modulus in the absence 
of pins in the weak field limit as noted by Fetter, Hohcn- 
berg and Pincus0,EEl. The purpose of this paper is 
to address the importance of this omission by investigat- 
ing the properties of the pinned and moving phases via 
3D molecular dynamics simulations of vortices interact- 
ing via a long-range potential in the presence of colum- 
nar defects. Our key result is that the softening of the 
shear modulus in the low field limit has important impli- 
cations for the phase diagram. We confirm the existence 
of a weakly pinned interstitial glass which becomes more 
weakly pinned as the field increases in agreement with 
previous studies of short-range interactions 4] . However 
we find a qualitatively new phase diagram for low fields 
where the interstitial glass melts near B = B^ and recrys- 
talizes at higher fields when the matching field is below 
a critical value. In this regime we find strong numeri- 
cal evidence for enhanced values of the critical current 
near three times the matching field. This questions the 
quantitative appropriateness of the mapping of vortex 
physics onto the physics of the superconductor-insulator 
transition without a properly detailed consideration of 
long-range interactions. 

We model the motion of vortices as coupled pancakes 
on neighboring, continuous planes separated by a dis- 



tance z, with the applied magnetic field aligned perpen- 
dicularly to the planes. The off-lattice simulation models 
the motion of vortices referenced by a 2D coordinate r 
under the influence of pinning and repulsive vortex inter- 
actions, and driving, line bending, and thermal forces: 



f T - 



0g({ ri ,r a ...}) 



(1) 



Here fjj = QqJ/c is the Lorentz force per unit length 
due to an applied current density J perpendicular to the 



magnetic field. r\i = 



cT-, 



is the Bardeen-Stephen vis- 



cous drag coefficient, with <£> = (hc/2e) the flux quan- 
tum and p n the normal state resistivity. The Langevin 
thermal force per unit length fa is normalized to set the 
rms vortex velocity via the equipartition theorem. The 
Hamiltonian H is the sum of the line tension for bend- 
ing, and vortex- vortex and vortex-disorder potentials per 
unit length constructed via London theory. The line ten- 
sion is given by e; = eln(K)eo with e = (m c /m a fc) 2 ,K = 
V£> e o = (^o/47rA) 2 , and A magnetic penetration depth. 
The vortex interaction is given by a sum of pairwise inter- 
actions V v - V (r) — eoK(r/X) with K a Hankel functionyj. 
We model the correlated extended defects as smooth 
parabolic traps of width R p and uniform depth V p = e;/4. 
Defects are randomly placed and aligned along the c-axis. 

For DC transport, in most cases the vortex mass per 
unit length mi is overall quite small in comparison with 
the other parameters in Eq. Q 18] and has thus usu- 
ally been neglected in p revious numerical studies in 3D 
of vortex dynamics 19]. However, in the case of su- 
perconductors near a Mott instability, such as the un- 
derdoped cuprates, the vortex mass can be immensely 
enhanced as the system approaches a superconductor- 
insulator transition 20] . We thus keep this term and 
simply relate the vortex mass to the mass of the highly 
renormalized electrons confined to the vortex core region, 
mi = m e f ftiiT^ 2 , with n the electron density and m e ft 
the in-plane effective mass of the electron renormalized 
by strong interactions. While we have found that the 
vortex mass has an impact on the magnitude of the crit- 
ical current, the relative J c for different vortex densities 
do not drastically depend on the choices made. 

Periodic boundary conditions are imposed in the 
planes to maintain constant global flux density, and open 
boundary conditions are employed along the c-axis. Tem- 
perature is chosen which is high enough to allow individ- 
ual vortices to be quickly accommodated to defects in 
the absence of a driving current but well below the glass 
temperature^. 

We measure all energies in units of the bare line bend- 
ing energy eo and measure all lengths in units of d = 4£. 
A natural time unit to is chosen to be irnr]id 2 /ei and 
the time step is further discretized in units of O.Olto for 
the simulations. The current is measured in units of the 
BCS depairing current Jq = $ock/12-\/3V 2 A 3 , and resis- 
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FIG. 1: Resistivity as a function of driving current for dif- 
ferent values of the applied magnetic field referenced to the 
matching field B$ for B^/H c2 x 10 3 = 7.32 (A) and 1.525 (B), 
respectively. The values of B/B<s> are indicated as follows: 
Panel (A): 0.3 (circles), 0.8 (squares), 1 (diamonds), 1.25 (up 
triangles), 1.48 (left triangles), and 2 (down triangles); and 
Panel (B): 0.5 (circles), 1 (squares), 1.5 (diamonds), 2 (up 
triangles), 3 (left triangles), and 4 (down triangles). Solid 
lines are guides. The insets show the low current resistivi- 
ties. Currents,resistivities are measured in units of Jo,Pbs 
respectively, as defined in the text. Color on-line. 



tivities in terms of the Bardeen-Stephen flux-flow resis- 
tivity pbs = B&o/ ! {c 2 rji). The other parameters used 
in the simulations are dictated by values appropriate 
for YBa 2 Cu307: m e ff = 5m e ,R p — 2£,£ = 17 A, z = 
12 A, k = 100, and e = 1/25, giving H c2 = 120T. 

Our simulations were performed using up to 40,000 
vortex segments in a 128c? x 128d square periodic cell 
containing 80 planes, where finite size effects were investi- 
gated and found to be minimal. Measurements are taken 
over an interval of 10 6 io after the system has reached a 
steady state after a typical time 3 x 10 5 t . We typically 
average the results over several hundred realizations of 
disorder particularly at smaller driving forces. We mea- 
sure the average vortex velocity in the direction of the 
Lorentz force corresponding to the voltage drop across 
the sample, and determine the resistivity p — £/ J. 

Our results for the resistivity p versus the applied cur- 
rent density J are shown for a series of vortex densities 
below and above the matching field for two different 
values of _B$ in Fig. (^V and B). All error bars are equal 
to the symbol size. For the pinned Bose-glass regime 
B < _B$ and for appreciable Z?$ (Fig. 1A), vortices are 
localized on separate columnar defects for small J un- 
til an abrupt transition to a moving regime ensues near 
J ~ O.Uo where the vortices simultaneously become un- 
pinned. For larger driving currents, vortices are in the 
flux flow regime and the resistivity approaches pss- The 
de-pinning transition occurs within a narrow range of 
currents corresponding to single vortex pinning. Conse- 
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FIG. 2: Normalized critical currents as a function of B/B$> 
for different densities of columnar defects. The density of 
defects is given by 10 3 x B$/H c2 = 1.525,1.906,3.05,4.88, 
and 7.32 for circles, squares, down triangles, diamonds, and 
up triangles, respectively. Color on-line. 



quently, as B increases, the depinning transition occurs 
for smaller J and significantly broadens as the effective 
role of disorder diminishes relative to the vortex lattice 
interaction, in agreement with previous simulations 
The resistivity increases monotonically with increasing 
B as shown in the inset of Fig. (1A). For values of B 
near and above the matching field, a low-current tail de- 
velops for small J and continues to increase as more and 
more interstitial vortices are not able to reside on colum- 
nar defects due to the repulsion from the pinned vortices. 
These interstitial vortices are caged by the pinned vor- 
tices forming an weakly pinned interstitial glass with a 
much reduced critical current J c , as suggested in Refs. 
0, 0- As i? is further increased, the interstitial glass 
melts into a interstitial liquid where channels of vortices 
flow around pinned vortices for arbitrarily small J and 
give a resistivity proportional to the fraction of intersti- 
tial to pinned vortices. We note that for our parameter 
choices we do not observe a Mott insulating phase as we 
always have interstitial vortices for B /_B$ = 1 due to the 
long-range vortex repulsion and high temperatures. 

Important differences however are observed for smaller 
values of £?$ (Fig. IB). For B < £>$ in this case, the tran- 
sition to depinning occurs at lower J and is substantially 
broadened compared to Fig. (1A). A non-monotonic de- 
pendence of the resistivity appears as shown in the in- 
set of Fig. (IB). For low driving forces the resistivity 
rises with increasing B for B/B§ < 1.5 and becomes 
greater than O.lpBS for B > S$ as the interstitial glass 
melts into an interstitial liquid as the density of vortices 
is increased. However at larger fields the resistivity de- 
creases and becomes less than Q.lpss for fields near three 
times the matching field. This indicates that for low B$, 
the Bose glass melts into an interstitial liquid near the 
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FIG. 3: Conjectured phase diagram for large (A), intermedi- 
ate (B), and small (C) matching field densities. 



matching field and recrystalizes into an interstitial glass 
at larger values of B/B$, in contradiction to the phase 
diagram proposed previously^]. 

We can make this more quantitative by defining a crit- 
ical current density J c as the value of the current density 
corresponding to a resistivity of 10% of pss- For a num- 
ber of different molecular dynamics runs of different B$, , 
the values we obtain for J c are shown in Fig. (2) as a 
function of B / B§ normalized to the values of J c deter- 
mined in the small B/B& limit. For large the critical 
current decreases monotonically with B / B<$, with a grad- 
ual fall-off near the matching field. For smaller values 
of B$> the fall-off near B^> becomes much more abrupt, 
suggesting that a magnitude of the kink in the melting 
curve would be dependent on the value of -B*, reconcil- 
ing previous experiments!! HHH[lJH[ll[ll[li. 
For even smaller B& , the fall-off is dramatic and the in- 
terstitial glass weakens appreciably and melts (J c = 0) 
near a critical value of B^/H c2 ~ 1.9 x 10 -3 . Remark- 
ably, the critical current resurrects for larger values of 
the magnetic field and has a broad peak near B = 3S$ 
for B<f,/H c2 x 10 3 = 1.906 before falling off again at still 
larger fields. This reproduces the peak near three times 
the matching field seen in experiments on Tl2Ba2Cu06+<5 
. However we see that still smaller values of B$ yields 
a critical current peak at larger values of BjBq, suggest- 
ing that this peak would be dependent on the number of 
defects. While we expect that the actual value of the crit- 
ical field B c might depend on our choice of defining J c , 
we do not expect the rentrant behavior into the intersti- 
tial glass to be qualitativelty changed. A proposed phase 
diagram which encompasses our results for different J5$ 



is given in Figure J3Jl. 

Using renormalization group arguments it was shown 
by Nelson and Seung that for clean systems the interac- 
tions between vortex lines are renormalized to zero near 
H c i and leads to the melting of the Abrikosov lattice as 
H — > H c i 01 • Thus we might expect that the melting of 
the interstitial glass for B — Bq, > + is similar to the 
melting of the vortex lattice near H c i 0, 0] . However 
for B < B§ vortex localization onto unoccupied colum- 
nar defects helps to reduce vortex fluctuations and one 
might expect the Bose glass to be stable near H c i as vor- 
tices occupy the strongest defects. Yet once B > £?$ the 
interstitial vortices feel a much reduced and screened de- 
fect interaction and suffer the downward renormalization 
of the shear modulus as in clean systems. This is borne 
out in our simulations for small defect concentrations. 
As the field increases our simulations indicate that the 
interstitial glass recrystalizcs. This can be viewed as a 
hardening of the shear modulus for increasing fields as 
the role of interactions increases. 

What leads to the recrystalization of the melted inter- 
stitials at larger values of B for small defect concentra- 
tions? At large vortex concentrations the Coulomb po- 
tential is effectively short-range due to screening and the 
effects of disorder arc also screened by flux line collisions, 
and indeed our results are consistent with prior simula- 
tions for contact repulsive potentials 0. The system fa- 
vors forming dislocations at larger fields and the domains 
are vortices are easily depinned. For smaller vortex con- 
centrations, new physics arises as screening is no longer 
effective and the bare long-range repulsive forces encour- 
age longer-range ordering concommitant with decreased 
prevalence of dislocation and subsequent enhanced vor- 
tex pinning. 

To check these ideas, we plot in Figure (@J time- 
averaged structure plots for B = 3B$ at two values 
of -B$. While long-range order is not seen even on the 
length scale of our simulations, it is clear that orienta- 
tional short-range order is more prevalent for the small 
-B$ results than the large ones, with subsidiary struc- 
ture peaks at the reciprocal lattice vectors up to 1/2 the 
height of the central (Q = 0) peak. When averaged over 
many disorder configurations, the orientational order is 
lost for both values of yet it is clear that the effective 
range of the vortex repulsion is longer for smaller _B$ . We 
would expect that for still smaller values of B$ positional 
order would grow for systems of fixed finite size L with 
the possible formation of Bragg-like peaks. Most likely 
the positional order would not be quantitatively relevant 
for L — > oo. 

In summary, we have presented numerical simulations 
of interacting vortex dynamics in the presence of colum- 
nar disorder as a function of Bj £?$ and Bq> j H C 2 . We find 
a monotonic decrease in the critical current for large val- 
ues of -B$,a with a sharp drop-off near the matching field, 
consistent with prior notions of the weakly pinned inter- 
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FIG. 4: Structure factor in the pinned phase at three times 
the matching field B / B& = 3 for 10 3 x B$/H c2 = 1.525 (left) 
and 7.32 (right) panels, respectively. Momentum is measured 
in units of 1/d. Color on-line. 



sitial glass emerging at higher fields. For smaller values 
of B$ however, we see abrupt melting at B = B& and 
recrystalization near B = 3B§, of the interstitial glass, 
suggesting modifications to the phase diagram due to the 
role of long-range interactions even for columnar disor- 
der. 
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